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The first experimental demonstration of an adaptive quantum state estimation (AQSE) is re- 
ported. The strong consistency and asymptotic efficiency of AQSE have been mathematically 
proven [J. Phys. A:Math. Gen. 39 12489 (2006)]. In this Letter, the angle of linear polariza- 
tion of single photons, or the phase parameter between the right and the left circularly polarization, 
is estimated using AQSE, and the strong consistency and asymptotic efficiency are experimentally 
verified. AQSE will provide a general useful method in both quantum information processing and 
metrology. 
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Quantum theory is inherently statistical. This entails 
repetition of experiments over a number of identically 
prepared quantum objects, for example, quantum states, 
if one wants to know the "true state" or the "true value" 
of the parameter that specifies the quantum state [1-4]. 
Such an estimation procedure is particularly important 
for quantum communication and quantum computation 
[5], and is also indispensable to quantum metrology [6- 
10]. In applications, one needs to design the estimation 
procedure in such a way that the estimated value of the 
parameter should be close to the true value (consistency), 
and that the uncertainty of the estimated value should be 
as small as possible (efficiency) for a given limited num- 
ber of samples. In order to realize these requirements, 
Nagaoka advocated an adaptive quantum state estima- 
tion (AQSE) procedure [11, 12], and recently Fujiwara 
proved the strong consistency and asymptotic efficiency 
for AQSE [13, 14]. 

In this letter, we report the first experimental demon- 
stration of AQSE using photons. The angle of a half 
wave plate (HWP) that initializes the linear polarization 
of input photons is estimated using AQSE. A sequence 
of AQSE is carried out with 300 input photons, and the 
sequence is repeated 500 times for four different settings 
of HWP. The statistical analysis of these results veri- 
fies the strong consistency and asymptotic efficiency of 
AQSE. Recently, it has been mathematically proven that 
the precision of AQSE outperforms the conventional state 
tomography [15]. It is thus expected that AQSE will pro- 
vide a useful methodology in the broad area of quantum 
information processing, communication, and metrology. 

Let us first explain AQSE in detail. For simplicity, 
we restrict ourselves to one-dimensional quantum statis- 
tical model S = {pg; 9 £ O (c R)}, a smooth parametric 
family of density operators on a Hilbert space ri hav- 
ing a one-dimensional parameter 0. Our aim is to es- 



timate the true value of 9 by means of a certain quan- 
tum estimation scheme. An estimator is represented by 
a pair (M, 6), where M = {M{x)\ x € X} is a positive 
operator-valued measure (POVM) that takes values on a 
set X, and 9 : X — > 9 is a map that gives the estimated 
value 9{x) from each observed data x £ X. The observed 
data x € X has probability density 



f(x;9,M) :=Tr p e M(x), 



(1) 



which depends on both the parameter 9 and the mea- 
surement M. 

In traditional statistics, it is often the case to con- 
fine our attention to unbiased estimators. An estimator 
(M, 9) is called unbiased if 



Eg [M, 9} 



(2) 



is satisfied for all 9 £ O, where Eg[ ■ ] denotes the expec- 
tation with respect to the density (1). It is well known 
[16] that an unbiased estimator (M, 9) satisfies the quan- 
tum Cramer- Rao inequality Vg[M, 9] > (Jg)~ , where 
Vg[ ■ ] denotes the variance, and Jg is the quantum Fisher 
information of the model iS defined by Jg := Tr pgLg, 
where Lg is the symmetric logarithmic derivative (SLD) 
defined by the self-adjoint operator satisfying the equa- 
tion ^ = 5 {Lgpg + pgLg) . 

In quantum statistics, however, it is regarded that un- 
biasedness is too restrictive a requirement, and we usu- 
ally weaken the condition to a "local" one. An estima- 
tor (M, 9) is called locally unbiased [17] at a given point 
9 € if the condition (2) is satisfied around 9 = 9 Q 
up to the first order of the Taylor expansion, that is, if 
Eg Q [M,9] = O and ^Eg[ M ^}\ e= g = 1 hold - Clearly, 
an estimator is unbiased if and only if it is locally un- 
biased at all 6 € O. A crucial observation is that an 
estimator (M, 9) that is locally unbiased at 9q also satis- 
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FIG. 1: Schematic of adaptive quantum state estimation. 
Photons are linearly polarized with a polarization direction 
determined by HWPO. The polarization is analyzed by HWP1 
and the polarizing beam splitter (PBS). The controller sets 
HWP1 to an angle calculated on the basis of the photon mea- 
surement results. 



lies the quantum Cramer-Rao inequality 

V 6o [M,6]> (JeJ" 1 (3) 

at 9 = #n, and that the lower bound in (3) is achiev- 
able for any one-dimensional quantum statistical model 
S. To put it differently, the best locally unbiased estima- 
tor (LUE) for the parameter 9 at 9 = 9q is the one that 
satisfies Vg [M, 6] = (Jg,,) -1 . 

Here we encounter a difficulty which often becomes 
the target of criticism: since the best LUE for estimating 
the parameter 9 depends, in general, on the unknown 
parameter 9 itself, the estimation strategy based on LUEs 
would be infcasiblc. In a different yet analogous context, 
Cochran [18] ingeniously described this kind of dilemma 
as follows: "You tell me the value of 9 and I promise to 
design the best experiment for estimating 9P 

To surmount this difficulty, Nagaoka [If, 12] advocated 
an adaptive quantum state estimation (AQSE) scheme 
as follows. Suppose that, by prior investigation of the 
quantum statistical model S, one has the list of optimal 
LUEs (M{ ■ ;O),0( ■ ;<?)) for each G 8. One begins 
with an arbitrary initial guess ^ £ 8, and applies the 
measurement M{ ■ ;#n) that is optimal at 9q. Suppose 
the data x\ is observed, one then applies the maximum 
likelihood method to the likelihood function L\(ff) = 
f(x\;9,M( ■ ;9q)), to obtain the next guess 9\. At 
stage n (> 2), one applies the measurement M( ■ ; 0,i-i), 
where # n _i is the maximum likelihood estimator (MLE) 
obtained at the previous stage. The likelihood function 
is then given by L n (9) := flLi f{xi]0,M( ■ 
where Xi is the observed data at stage i, and one ob- 
tains the nth MLE 9 n that maximizes L n {9). It is quite 
natural to expect that the sequence 9 n of MLEs would 
converge to the true value of the parameter 9. In fact, 
under certain regularity conditions, it can be shown that 



the sequence 9 n is strongly consistent and asymptotically 
efficient [13, 14]. 

Now let us discuss the implementation of AQSE using 
photons (Fig. 1). Here the unknown parameter is the 
angle 9 of HWPO, which determines the phase </> between 
right and left circularly polarizations of input photons by 
the relation <j) = A9. An arbitrary linear polarization can 
be described using right and left circular polarizations as 
follows: 

IV.) = ±(\R) +e i * \L)) = cos(^) \H) + S in(^) \V) . (4) 

By changing the angle of the half wave plate (HWP1), 
we can adjust the measurement basis. For such measure- 
ment, the POVM having optimal estimation capability 
is given by 

M{9) = (M(l; 0), M(2; 9)) = (|0 (£| , / - |£) (£|) , (5) 

where (f | = (cos (29 + f ) , sin (20 + f )) . By applying 
the POVM M(9) to the input state \ip(0)) := \ one 
obtains the probability distribution on X := {1, 2} which 
is isomorphic to the fair coin flipping. 

The drawback to realizing this measurement is that 
the optimal POVM M(ff) depends on the unknown value 
of the parameter 9 [19]. We can avoid this drawback by 
adopting an AQSE as follows. We begin by setting the 
initial log-likelihood function to be Iq(8) = 0, and then 
start inputting and detecting photons one by one. For 
?ith photon, we apply the measurement M(6 n —i) which 
depends on the latest MLE 0„_i. Let x n E X be the 
outcome indicating which detector has been lit. The log- 
likelihood function is then updated by the formula 

l n (0) ~ ^(6) + log (tP(6)\M(x n ; 6^)^(6)), (6) 

and the nth MLE is given by 9 n = arg maxe l n (9) . Let 
us denote the true value of the parameter 9 by 9 l . It is 
known [13, 14] that the sequence 9 n of MLEs converges 
to the true value 9 l with probability one (strong con- 
sistency) and that the distributions of the random vari- 
ables \pn(9 n — 9 r ) converge to the normal distribution 
N(0, Jgt 1 ) (asymptotic efficiency), where Jg denotes the 
quantum Fisher information of the parameter 0, which 
turns out to be 16 for our model (4). 

The experimental setup is shown in Fig. 2(a). Single 
photons at 780nm are generated from a heralded sin- 
gle photon source [20], consisting of a CW diode pump 
laser (wavelength: 402 nm) and a 3 mm long BBO crys- 
tal (Type I). A pair of a signal photon (780 nm) and 
a trigger photon (830 nm) is created via spontaneous 
parametric down conversion. The detector (DT, SPCM- 
AQR, Perkin Elmer) after an interference filter (IF1, cen- 
ter wavelength 830nm) outputs an electric pulse (width 
30ns) when it detects a trigger photon and the electric 
pulse heralds the generation of a signal photon, which 
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FIG. 2: (a) Schematic of the experimental setup. (b)(c) An 
example showing the update of a log-likelihood function. The 
second term log {ip(9)\ M(x n ',O n -i) \ip{0)) in eq. (6) is shown 
in panel (b), and the updated l„(0) is shown in panel (c). The 
blue arrows indicate the true value . 



FIG. 3: (a) Trajectories of estimated HWPO angles against 
the number n of photons for r — 500 repetitions is shown in a 
three dimensional plot, (b) The first 10 curves are superposed 
in a two dimensional graph. 



is coupled to a polarization maintaining fiber (PMF) af- 
ter an interference filter (IF2, center wavelength 780 nm, 
width 4 nm). The polarization of photons are then initial- 
ized to be horizontal using a polarizer (extinction ratio 
10 -5 ). The target parameter 0* was set using HWPO. 
The polarization state of the photon was analyzed by 
HWP1 and a polarizing beam splitter (PBS). After pass- 
ing through the PBS, photons are guided to single pho- 
ton detectors (DO and Dl, SPCM-AQR, Pcrkin Elmer) 
on each PBS output port. The outputs of single photon 
detectors are gated by the rise of the heralding signal 
and connected to the "first-come discriminator," consist- 
ing of a home-made electric circuit. When the discrim- 
inator receives the first signal from one of the detectors 
(DO or Dl) after the measurement for (n — l)th photon 
starts, the discriminator informs which detector has been 
clicked. The minimum pulse interval of 2.5ns can be dis- 
criminated. Note that the discriminator ignores the case 
when it receives the pulses from both the detectors within 
2.5ns. The angle of HWP1 for measuring the nth photon 
is determined by calculating the discretized MLE 9 n , the 
maximizer of the log-likelihood function (6) chosen from 
among the 10000 points that divide the domain [0, 7r/2) 
of the parameter 9 into equal parts (Figs. 2(b) and 2(c)). 
When the change of HWP1 angle is completed, the mea- 
surement for the next (nth) photon will be started. In 
a sequence of AQSE, the above mentioned procedure is 
carried out up to 300 input photons (n=300). For four 
different HWPO angles 9 = 0,30,60, and 78.3 [deg], we 
repeated the sequence for 500 times (r=500). 



Let us first observe the strong consistency for the se- 
quence 9 n of MLEs for the parameter 9 of HWPO. Fig. 
3 (a) shows 500 trajectories of estimated HWPO angle 
9 n against the number n of photons when the true value 
0* of the parameter is set to be 60 degree. The curves 
correspond to independent runs of adaptive estimation. 
Evidently, each curve of 9 n approaches the true value 
0', which is in accord with the mathematical result that 
9 n — > 9 l almost surely as n — > oo, even though the curves 
are dissimilar to each other reflecting the genuine statis- 
tical nature of quantum system. The convergence to the 
true value is clear in Fig. 3(b) where first 10 trajectories 
in Fig. 3(a) are superposed. 

We next test the hypothesis that the MLE 9 n follows a 
normal distribution for large n. More concretely, we will 
investigate if the random variable y/nJg {9 n — 9) follows 
the standard normal distribution N(0, 1), i.e., %/nJg (9 n — 
9) ~ N(0, 1), where 9 is the sample average of MLEs 9 n 
over sufficiently many independent trials. A goodness of 
fit test [21] was carried out as follows: 

1) The real axis was divided into 23 intervals (bins) 
{Ib}bLo, where Ji,...,/2i are disjoint partitions of the 
interval [—3.5, 3.5] of equal width, and Iq = (— oo, —3.5), 
I22 = (3.5, +00). In reality, these bins were slightly 
shifted by 5/10000, where 5 := % AOe7r/20000 is the 
scaled resolution of the estimator 0„, so that the data 
-y/nJe (0„ — 9) did not fall on the boundaries of the bins. 

2) The test-statistic X 2 := Y%L ^"rp^ was calcu " 
lated, where Nb is the number of observed data which 
fell into 6th bin, pb the theoretical probability of falling 
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FIG. 4: Histogram of the observed data obtained by r = 500 
independent experiments of adaptive estimation scheme, each 
using n = 300 photons. These histograms were taken for four 
different true values of (a) [deg], (b) 30 [deg], (c) 60 [deg] 
and (d) 78.3 [deg]. 



a datum into 6th bin under the null hypothesis N(0, 1), 
and r the number of repetitions of adaptive estimation 
procedure. 

3) The test-statistic X 2 was analyzed using the chi- 
square distribution x| 3 _ p of degree 23 — p, where p = 2 
degrees of freedom ought to be subtracted because of the 
normalization and the use of sample average 9. 

Figure 4 shows the histogram of the observed data ob- 
tained by r = 500 independent experiments of adaptive 
estimation scheme, each using n = 300 photons. The 
true values 6>* of the parameter 9 of HWP0 are set to 
be 0, 30, 60, and 78.3 degrees. The density function of 
the standard normal distribution N(0, 1) is also plotted 
as the solid curve. All the experimental data agree with 
the standard normal distribution. To be precise, the val- 
ues of the test statistic X 2 arc (a) 16.8 (b) 15.7 (c) 12.8 
(d) 16.2, and the null hypothesis is accepted with 10% 
significance level in each case. 

Having obtained the strong evidence that the distri- 
bution of the MLE has converged quite well to a nor- 
mal distribution at n = 300, we finally proceed to the 
estimation of confidence intervals [21] for the mean p 
and variance v, assuming that \pn(9 ri — p) ~ N(0,v). 
The confidence intervals for p and v arc obtained by 
the standard procedure based on the statistical laws that 
(0 - p) ~ T r _r and ^ V ~ x 2 r -v Here V is the 

unbiased variance of MLEs 9 n over r trials, and T r —\ the 
^-distribution of degree r — 1. 

Table 1 summarizes the results for r = 500 with 90% 
confidence level. Recall that the asymptotic efficiency 
asserts that p ~ 9 t and v ~ Jgt (= 0.0625). Since the 
precision of the present experiment is about ±0.2 degree 
[22], we conclude that the estimated values of p and v 
listed in Table I are in excellent agreement with the the- 



TABLE I: Confidence intervals for the mean p and 
the variance v. CL means confidence level. 



l [deg] 


p [deg] (90% CL) 


v (90% CL) 


0.0 


-0.15 ± 0.06 


[0.054, 0.067] 


30.0 


29.90 ± 0.06 


[0.055, 0.067] 


60.0 


60.00 ± 0.06 


[0.056, 0.068] 


78.3 


78.27 ± 0.06 


[0.055, 0.068] 



orctical values. 

It should be noted that the purpose of our AQSE is 
completely different from 'adaptive measurements' pro- 
posed by Berry and Wiseman [23]. Their scheme was 
devised to estimate the phase difference between the two 
arms of an interferometer using a special ./V-photon two- 
mode state, approximating the canonical measurement 
proposed by Sanders and Milburn [24], and is not ap- 
plicable to general quantum state estimation problems. 
By contrast, our AQSE is a general-purpose estimation 
scheme applicable to any quantum statistical model using 
n identical copies of an unknown state. AQSE may also 
be used in verifying the achievability of the Cramer- Rao 
version of the Heisenberg limit 0(1/N 2 ) [25] by apply- 
ing the scheme to the n-i.i.d. extension pf n of an TV- 
photon phase-shift model pe on T-L ~ (C 2 )® N . (See also 
[26] for estimating a unitary channel under noise.) Inci- 
dentally, AQSE is based on the Cramer-Rao type point 
estimation theory and is free from the choice of a priori 
distribution which matters in Bayesian statistics such as 
adaptive Bayesian quantum tomography [27]. 

In summary, we have verified both the strong con- 
sistency and asymptotic efficiency of AQSE by experi- 
mentally estimating the angle of linear polarization of 
photons. Since AQSE has been mathematically proven 
to outperform the conventional estimation scheme such 
as the state tomography [15], we plan to apply AQSE 
to multi-parameter cases and compare the performance 
with other protocols using fixed measurement basis [28] . 
It will also be intriguing to apply AQSE to enhance the 
performance of quantum metrological experiments beat- 
ing the standard quantum limit [6-9] . 
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